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Outline: 

• What is the pion distribution amplitude (p n (x)7 

• Nonperturbative part: How to obtain ^^{x) from QCD sum rules; 

• Perturbative part: NLO light-cone sum rules =>- CLEO experiment on F 11 V {Q 2 ) =>• 
constraints on <f n (x) and A 2 ; 

• Perturbative addition: Diffractive dijet production (E791 data); 

• Perturbative addition: Pion electromagnetic form factor (CEBAF data); 

• Conclusions. 

The main object of this talk is the pion distribution amplitude (DA), which 
can be defined through the matrix element of a nonlocal axial current on the light cone 



(0\d(z)^ 5 E(z,0)u(0)\n(P)) 



iUP„ / dxe^ <^(x )M 2 



MP) 



which is explicitly gauge-invariant due to the presence of the Fock-Schwinger connector 
E(z,0) = p e i9 ^ A ^ dT \ The physical meaning of this object is quite transparent: 

It is the amplitude for the transition of 
the physical pion ir(P) to a pair of va- 
lence quarks u and d, separated at light- 
cone (see graphical image to the right), 
with momentum fractions xP and xP, 
correspondingly (here x = 1 — x) . 

This object inevitably appears in applying perturbative QCD to hard processes with pi- 
ons in the initial or the final state as a result of QCD factorization theorems E3 
and it includes nonperturbative information about the physical pion. It has the following 
properties: 




normalized to unity f Q dx(p n (x, // 2 ) = 1; 

x x symmetric: ip n (x,fi 2 ) = <f n (x, /z 2 ); 

obeys the ER-BL evolution equation with respect to /i 2 ; 

in the 1-loop approximation ^(x; /i 2 — > oo) = ip as (x) = 6x(l — x). 



1 Talk presented at the 13th International Seminar on High Energy Physics "Quarks-2004" , Pushkinskie 
Gory, Russia, May 24-30, 2004. It is based on results obtained in collaboration with S. Mikhailov, 
K. Passek-Kumericki, W. Schroers, N. G. Stefanis 
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It is convenient to represent the pion DA as an expansion in terms of Gegenbauer 
polynomials C^f 2 (2x — 1), being the 1-loop eigenfunctions of the ER-BL kernel: 



^(x;/i 2 ) = ^ as (x)[l+a 2 (^ 2 )C 2 3/2 (2x-l)+a 4 (/i 2 )Q 3/2 (2x-l) + ... 



(2) 



That means to transfer all the /^-dependence of the pion DA into the Gegenbauer coef- 
ficients {a2(/i 2 ), a^(fi 2 ), . . .}. This scheme can be effectively applied at the 2- loop level as 
well HE]. 
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Figure 1: (a) The bunch of pion DAs extracted from NLC QCD sum rules. Parameters of the bold- 
faced curve are 03 = +0.188 and = —0.130. (b) The results for the inverse moment (a:" 1 )^ as a 
function of the Borel parameter M 2 obtained using a special model-independent sum rule. The shaded 
area corresponds to the 10%- variation of the threshold parameter so- Dashed straight lines show the 
allowed window for (x" 1 )^. 

In order to obtain the pion DA in the theory, one is obliged to use some 
nonperturbative approach. Historically, the first nontrivial model has been constructed 
by Chernyak and Zhitnitsky (CZ) [H] using the standard QCD sum rule approach and 
estimating the first two moments of the pion DA: (£ 2 ),r and {£, 4 ) n - After that, Mikhailov 
and Radyushkin realized that in doing so CZ highly overestimated these moments and 
suggested to use instead the non-local condensate (NLC) approach [7] . We have used the 
NLC QCD sum rules and obtained the first five moments of the pion DA, (C, 2N )tt with 
N — 1, ...,5. Just for illustration, we present here the simplest scalar condensate of the 
used NLC model: 

(q(0)q(z)) = (g(0)g(0)} e~«/ 8 . (3) 

This model is determined by a single scale parameter A 2 = (k 2 ) characterizing the average 
momentum of quarks in the QCD vacuum. It has been estimated in QCD SRs and on 
the lattice: 

r 0.4 ± 0.1 GeV 2 [ QCD SRs ] 
A 2 = < 0.5 ± 0.05 GeV 2 [ QCD SRs [H] ] (4) 
[ « 0.4 - 0.5 GeV 2 [ Lattice [TUl EE] ] 

NLC sum rules for the pion DA produce |I2| a "bunch" of self-consistent 2-parameter 
models at /a 2 ~ 1 GeV 2 : 

(p n (x) = </? as (x) [l + a 2 C 2 3/2 (2x-l)+a 4 C 4 3/2 (2x-l)l . (5) 
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Figure 2: (a) The bunches of pion DAs extracted from NLC QCD sum rules in the (02, a4)-plane for three 
values of the nonlocality parameter A 2 , (b) Comparison of the NLC-bunch evolved to /1 2 = 5.76 GeV 2 
with the CLEO data constraints for A 2 = 0.4 GeV 2 . The la- and the 2<r-contours are shown in dashed 
and solid lines, correspondingly. 



For the most favorite value of the vacuum nonlocality parameter A 2 = 0.4 GeV 2 we have 
the bunch of pion DAs presented in Fig.^. By self-consistency we mean that the value 
of the inverse moment for the whole bunch (^ x - l ^ nnch = 3.17±0.10 is in agreement with 
the independent estimation from the special sum rule, (x _1 )^ R = 3.30 ± 0.30, see Fig.[T)3. 

We also extract the corresponding bunches for two other values of A 2 , = 0.5 GeV 2 and 
A 2 = 0.6 GeV 2 , and show the results as allowed regions in the (02, a^-plane in Fig.|2k- 

NLO light-cone sum rules (LCSR) and the CLEO data on F 77 * 7T (Q 2 ) allow one 
to obtain constraints on f n (x) directly from the experimental data. A natural question 
arises: Why does one need to use LCSRs? The answer is that for Q 2 ^> m 2 , q 2 <C m 2 
pQCD factorization is valid only in leading twist but higher twists are also of impor- 
tance |13j . The reason is quite evident: if q 2 —>■ one needs to take into account the 
interaction of a real photon at long distances of order of 0(l/y / g 2 ). To account for long- 
distance effects in perturbative QCD, one needs to introduce a light-cone DA of a real 
photon. In the absence of reliable information about the photon DA, Khodjamirian [Tlj 
suggested to use the LCSR approach, which effectively accounts for long-distances effects 
of a real photon, using the quark-hadron duality in the vector channel and a dispersion 
relation in q 2 : 
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with So — 1.5 GeV 2 - effective threshold in vector channel, M 2 - Borel parameter (0.5 — 0.9 
GeV 2 ). We revised the NLO LCSR approach of ^3] in performing the CLEO data analysis 
along the following lines |16j : 

• An accurate NLO evolution for both ip(x,Ql xp ) and a s (Q 2 xp ), taking into account 
heavy quark thresholds. 

• The relation between the "nonlocality" scale and the the twist-4 magnitude #t w _ 4 ~ 
A 2 /2 was used to re-estimate <^ w _ 4 = 0.19 ± 0.02 at A 2 = 0.4 GeV 2 . 



Constraints on (a; 1 ) 7r from the CLEO data. 
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Figure 3: (a) Comparison of the NLC-bunch evolved to /i 2 = 5.76 GeV 2 with the CLEO data constraints 
for A 2 = 0.4 GeV 2 . The la-, 2a- and 3<r-contours are shown as dashed, solid and dash-dotted contours. 
For details see in the text, (b) Comparison of BMS, CZ and asymptotic DAs at the QCD sum-rule scale 
fj, 2 ~ 1 GeV 2 with the CLEO data constraints for A 2 = 0.4 GeV 2 in terms of rotated axes ai — 0,4 and 
a 2 + 04. The la- and 2cr-contours are shown as dashed and solid lines. 



As a result, we have obtained reasonable agreement of our bunch with the CLEO data 
for A 2 = 0.4 GeV 2 , see Fig.[2b (with (♦) = asymptotic DA, (X) = BMS model, (■) = CZ 
DA, and (+) corresponds to the best-fit point). 

In order to make our conclusions more valuable, we have adopted a 20% uncertainty 
in the magnitude of the twist-4 contribution, 5f w -4 = 0.19 ± 0.04 GeV 2 , and produced 
new la-, 2a- and 3<r-contours dictated by the CLEO data [T7] . see Fig.Ot in parallel with 
available 2-Gegenbauer models: asymptotic DA, BMS model, CZ DA (they are shown in 
the same manner as in Fig.|2b); three instanton-based models, viz., (it) [TS], ▲ |19j . and 
(♦) (using in this latter case m q = 325 MeV, n — 2, and A = 1 GeV) [20], and a recent 
transverse lattice result (▼) [21] ■ We see that even with a 20% uncertainty in the twist-4, 
the CZ DA is excluded at least at the 4<r-level, whereas the asymptotic DA - at the 
3<r-level. Our bunch is mainly inside the la-region and other nonperturbative models are 
near the 3<r-boundary. 

We also plot the CLEO data in the plane (X, Y) with X = ci2 — 04 and Y = a<i + = 
(3 ;_1 ) 7 r/3 — 1, where the Gegenbauer coefficients a 2 and 0,4 refer to the NLC sum-rule 
scale /i 2 ~ 1 GeV 2 . The result is shown in Fig.|Hb, where the comparison of the CLEO 
data constraints directly with the model-independent bound i(x _1 )| R — 1 = 0.1 ± 0.1 
from the NLC QCD sum rule (shaded strip in figure) is done. Again we see a good 
agreement of a theoretical "tool" of different origin with the CLEO data. Here, we should 
also mention other estimations of the pion DA inverse moment. Bijnens&Khodjamirian 
produced an estimate KoT 1 )^ - 1 = 0.24 ± 0.16 using data on the pion electromagnetic 
form factor in the LCSR approach [22], whereas Ruiz Arriola&Broniowski obtained in 
their model of the pion DA with an infinite number of Gegenbauer harmonics the result 
^(ar 1 )^ - 1 = 0.25 ±0.1 [2S[. 

To finish our discussion about the CLEO data constraints in the NLO LCSR approach, 
we show in Fig.E^i the plot of Q 2 F 1 * 1 ^ 1T (Q 2 ) for our bunch (shaded strip), CZ DA (upper 
dashed line), asymptotic DA (lower dashed line), and two instanton-based models (dot- 
ted [IB] and dash-dotted [21] lines) in comparison with the CELLO and the CLEO data. 
We see that the BMS bunch describes rather well all data for Q 2 > 1.5 GeV 2 . 

Diffractive Dijet Production: What can add the E791 data to our analysis? The 
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Figure 4: (a) 7*7 — > 7r Transition form factor in comparison with the CELLO (♦) 24 and the CLEO 
(a) data. For details see in the text, (b) Comparison of the asymptotic DA (solid line), CZ DA 
(dashed line), and the BMS "bunch" of pion DAs (strip) with the E791 data (♦) 

diffractive dijet production in 7r + A collisions has been suggested as a tool to extract 
the profile of the pion DA by Frankfurt et al. in 1993 28J. They argued that the jet 
distribution with respect to the longitudinal momentum fraction has to follow the quark 
momentum distribution in the pion and hence provides a direct measurement of the pion 
DA. As it was shown just recently in |2H] (see also (201), beyond the leading logarithms 
in energy this proportionality does not hold. Braun et al. found that the longitudinal 
momentum fraction distribution of the jets for the non-factorizable contribution turns out 
to be the same as for the factorizable contribution with the asymptotic pion distribution 
amplitude. We have used this convolution approach of Braun et al. to estimate the 
distribution of jets in this experiment for our bunch of pion DAs in comparison with the 
asymptotic and the CZ DAs [T7j. Results are shown in Fig.HjD. It is interesting to note 
that the corresponding \ 2 values 9X6 '. clS 12.56; CZ - 14.15; BMS - 10.96 (accounting 
for 18 data points). The main conclusion from this comparison: all three DAs are 
compatible with the E791 data. Hence, this experiment cannot serve as a safe profile 
indicator. 
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Figure 5: Comparison of the asymptotic (dotted line), the CZ (dashed line), and the BMS DAs (solid 
line) . 

Let us say a few words about similarities and differences between the CZ and BMS 
DAs. Both are two-humped, but the CZ DA is strongly end-point enhanced, whereas the 
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BMS DA is end-point suppressed! And the reason for this behaviour is physically evident: 
nonlocal quark condensate reduces pion DA in the small x region and enhances in the 
vicinity of the point x ~ 0.2. In order to keep the norm equal to unity, it is forced to have 
in the central region some reduction as well. 

Pion electromagnetic form factor: How well is the BMS bunch in comparison 
with the JLab data on the pion form factor? We have calculated the pion form factor in 
analytic NLO pQCD 



pFact-WI 



(Q 2 ;/4) 



(6) 



with taking into account the soft part F^ D (Q 2 ) via the local duality approach and the 
factorized contribution 



^Fact-WI 



(Q 2 ;/4) 



Q 2 



2s 2-loop + Q 2 



(7) 



has been corrected via a power-behaved pre-factor (with Sq" 1o ° p « 0.6 GeV 2 ) in order to 
respect the Ward identity at Q 2 = and preserve its high-Q 2 asymptotics. In our analysis 
Fj act (Q 2 ; has been computed to NLO [USES], using Analytic Perturbation Theory [3U 
EH] and trading the running coupling and its powers for analytic expressions in a non- 
power series expansion, i.e., 



K act (Q 2 ;/4)] 



MaxAn 



71 



with and 2 (a*r) being the 2-loop analytic images of a s (Q 2 ) and (a s (Q 2 )) 2 , corre- 
spondingly (see [3T] for more details), whereas JF^ (Q 2 ) and T^ h °{Q 2 ] are the LO 
and NLO parts of the factorized form factor, respectively. The result of this analysis is 
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Figure 6: Pion electromagnetic form factor in comparison with the JLab (♦) and Bebek et al. 
(a) data, (a) The solid line represents results obtained with /x^ = 1 GeV 2 , the dashed line - with 
/ip = Q 2 , the dotted line - with the BLM scale, and the dash-dotted line — in the ay-scheme, (b) 
Predictions based on the BMS "bunch" of pion DAs (strip) and the asymptotic DA (dashed lines). The 
green strip contains the NLC QCD sum-rule uncertainties (due to the BMS bunch) and scale-setting 
ambiguities at the NLO level (in the case of the asymptotic DA these ambiguities are represented by two 
dashed lines). 



presented in Fig. El where we show F n (Q 2 ) for the BMS "bunch" and using the "Maximally 
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Analytic" procedure, which improves the previously introduced "Naive Analytic" one. 
The new procedure with the analytic running coupling and analytic versions of its pow- 
ers gives us practical independence of the scheme/scale setting (see Fig.Ek and the figure 
caption for details) and provides results in a rather good agreement with the experimental 
data jHHIEZI- We see that the form-factor predictions are only slightly larger than those 
resulting when using the asymptotic DA (see Fig.EJa). 

Conclusions. 

• The QCD sum rule method with NLCs for the pion DA gives us admissible sets 
(bunches) of DAs for each X q value. 

• The NLO LCSR method produces new constraints on the pion DA parameters 
(02,04) in conjunction with the CLEO data. 

• Comparing NLC sum-rule results with the new CLEO constraints allows us to fix 
the value of QCD vacuum nonlocality \ 2 q = 0.4 GeV 2 . 

• The corresponding bunch of pion DAs agrees well with the E791 data on diffractive 
dijet production and with the JLab F(pi) data on the pion electromagnetic form 
factor. 

• Analytic perturbation theory with non-power NLO for the pion form factor dimin- 
ishes scale-setting ambiguities already at the NLO level, rendering still higher-order 
corrections virtually superfluous. 
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